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Abstract. 

The Strong CP Problem and its resolution through the existence of an axion are briefly reviewed. The combined constraints 
from accelerator searches, the evolution of red giants and the duration of the SN 1987a neutrino pulse require the axion mass 
m a to be less than 3 • 1CU 3 eV. On the other hand, the requirement that the axion does not overdose the universe implies a 
lower bound on m a of order 1CU 6 eV. This lower bound can, however, be relaxed in a number of ways. If m a is near the lower 
bound, axions are an important contribution to the energy density of the universe in the form of (very) cold dark matter. Dark 
matter axions can be searched for on Earth by stimulating their conversion to microwave photons in an electromagnetic cavity 
permeated by a magnetic field. Using this technique, limits on the local halo density have been placed by the Axion Dark 
Matter experiment. I'll give a status report on ADMX and its upgrade presently under construction. I'll also report on recent 
results from solar axion searches and laser experiments. 
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INTRODUCTION 

The standard model of elementary particles has been 
an extraordinary breakthrough in our description of the 
physical world. It agrees with experiment and observa- 
tion disconcertingly well and surely provides the founda- 
tion for all further progress in our field. It does however 
present us with a puzzle. 

Indeed the action density includes, in general, a term 
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where G^ v are the QCD field strengths, g is the QCD 
coupling constant and 9 is a parameter. The dots rep- 
resent all the other terms in the action density, i.e. the 
terms that lead to the numerous successes of the stan- 
dard model. Eq. (1) perversely shows the one term that 
isn't a success. Using the statement of the chiral anomaly 
[1], one can show three things about that term. First, that 
QCD physics depends on the value of the parameter 9 
because in the absence of such dependence QCD would 
have a Ua ( 1 ) symmetry in the chiral limit, and we know 
QCD has no such Ua(1) symmetry [2]. Second, that 9 
is cyclic, that is to say that physics at 9 is indistinguish- 
able from physics at 9 + In. Third, that an overall phase 
in the quark mass matrix m q can be removed by a redef- 
inition of the quark fields only if, at the same time, 9 
is shifted to 9 — argdetm^. The combination of standard 
model parameters 9 = 9 — argdetm^ is independent of 
quark field redefinitions. Physics therefore depends on 9 
solely through 9. 



Since physics depends on 9, the value of 9 is deter- 
mined by experiment. The term shown in Eq. (1) violates 
P and CP. This source of P and CP violation is incompat- 
ible with the experimental upper bound on the neutron 
electic dipole moment [3] unless |0| < 10~ 10 . The puz- 
zle aforementioned is why the value of 9 is so small. It is 
usually referred to as the "Strong CP Problem". If there 
were only strong interactions, a zero value of 9 could 
simply be a consequence of P and CP conservation. That 
would not be much of a puzzle. But there are also weak 
interactions and they, and therefore the standard model as 
a whole, violate P and CP. So these symmetries can not 
be invoked to set 0=0. More pointedly, P and CP vio- 
lation are introduced in the standard model by letting the 
elements of the quark mass matrix m q be arbitrary com- 
plex numbers [4]. In that case, one expects argdetm^, and 
hence 9, to be a random angle. 

The puzzle is removed if the action density is instead 
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where a(x) is a new scalar field, and the dots repre- 
sent the other terms of the standard model. f a is a con- 
stant with dimension of energy. The aG ■ G interaction in 
Eq. (2) is not renormalizable. However, there is a recipe 
for constructing renormalizable theories whose low en- 
ergy effective action density is of the form of Eq. (2). 
The recipe is as follows: construct the theory in such a 
way that it has a (7(1) symmetry which (1) is a global 



symmetry of the classical action density, (2) is broken by 
the color anomaly, and (3) is spontaneously broken. Such 
a symmetry is called Peccei-Quinn symmetry after its in- 
ventors [5]. Weinberg and Wilczek [6] pointed out that a 
theory with a £/pq(1) symmetry has a light pseudo-scalar 
particle, called the axion. The axion field is a(x). f a is 
of order the expectation value that breaks f/pQ(l), and is 
called the "axion decay constant". 

In the theory defined by Eq. (2), 9 — ^ — detargm,, 
depends on the expectation value of a(x). That expecta- 
tion value minimizes the effective potential. The strong 
CP problem is then solved because the minimum of the 
QCD effective potential V(B) occurs at 9 = [7]. The 
weak interactions induce a small value for 9, of order 
10~ 17 [8], but this is consistent with experiment. 

The notion of Peccei-Quinn (PQ) symmetry may seem 
contrived. Why should there be a 1/(1) symmetry which 
is broken at the quantum level but which is exact at the 
classical level? One should keep in mind, however, that 
the reasons for PQ symmetry may be deeper than we 
know at present. String theory contains many examples 
of symmetries which are exact classically but which 
are broken by anomalies, including PQ symmetry. Also, 
within field theory, there are examples of theories with 
automatic PQ symmetry, i.e. where PQ symmetry is a 
consequence of just the particle content of the theory 
without adjustment of parameters to special values. 

The first axion models had f a of order the weak inter- 
action scale and it was thought that this was an unavoid- 
able property of axion models. However, our conference 
chairman J.E. Kim and others pointed out [9, 10] that the 
value of f a is really arbitrary, that it is possible to con- 
struct axion models with any value of f a . J.E. Kim also 
pointed out [11] that a value of f a far from any previ- 
ously known scale need not lead to a hierarchy problem 
because PQ symmetry can be broken by the condensates 
of a new technicolor-like interaction. 

The properties of the axion can be derived using the 
methods of current algebra [12]. The axion mass is given 
in terms of f a by 
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All the axion couplings are inversely proportional to f a . 
For example, the axion coupling to two photons is: 
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(4) 



Here E and B are the electric and magnetic fields, a is 
the fine structure constant, and gy is a model-dependent 
coefficient of order one. gy = 0.36 in the DFSZ model 
[10] whereas g r = -0.97 in the KSVZ model [9]. The 
coupling of the axion to a spin 1/2 fermion / has the 
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FIGURE 1. Ranges of axion mass m a , or equivalently axion 
decay constant f a , which have been ruled out by accelerator 
searches, the evolution of red giants, the supernova SN1987a, 
and finally the axion cosmological energy density. 



form: 
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where gf is a model-dependent coefficient of order 
one. In the KSVZ model the coupling to electrons is 
zero at tree level. Models with this property are called 
'hadronic'. 

The axion has been searched for in many places but 
not found [13]. Fig. 1 summarizes the constraints. Axion 
masses larger than about 50 keV are ruled out by par- 
ticle physics experiments (beam dumps and rare decays) 
and nuclear physics experiments. The next range of axion 
masses, in decreasing order, is ruled out by stellar evo- 
lution arguments. The longevity of red giants rules out 
200 keV >m a > 0.5 eV [14, 15] in the case of hadronic 
axions, and 200 keV > m a > 10~ 2 eV [16] in the case 
of axions with a large coupling to electrons [g e = 0(1) in 
Eq. 5] . The duration of the neutrino pulse from supernova 
1987a rules out 2 eV > m a > 3 • 10~ 3 eV [17]. Finally, 
there is a lower limit, m a > 10~ 6 eV, from cosmology 
which will be discussed in detail in the next section. This 
leaves open an "axion window": 3 • 10~ 3 > m a > 10~ 6 
eV. Note however that the lower edge of this window 
(10~ 6 eV) is much softer than its upper edge. 



AXION COSMOLOGY 

The implications of the existence of an axion for the his- 
tory of the early universe may be briefly described as 
follows. At a temperature of order f a , a phase transition 
occurs in which the C/pg ( 1 ) symmetry becomes sponta- 
neously broken. This is called the PQ phase transition. 
At these temperatures, the non-perturbative QCD effects 
which produce the effective potential V(d) are negligi- 
ble [18], the axion is massless and all values of (a(x)) 
are equally likely. Axion strings appear as topological 
defects. One must distinguish two scenarios, depending 
on wether inflation occurs with reheat temperature lower 
(case 1) or higher (case 2) than the PQ transition tem- 
perature. In case 1 the axion field gets homogenized by 
inflation and the axion strings are 'blown' away. 

When the temperature approaches the QCD scale, the 
potential V(9) turns on and the axion acquires mass. 
There is a critical time, defined by m a (ti)t\ = 1, when 
the axion field starts to oscillate in response to the turn- 
on of the axion mass. The corresponding temperature 
T\ ~ 1 GeV [19]. The initial amplitude of this oscilla- 
tion is how far from zero the axion field lies when the 
axion mass turns on. The axion field oscillations do not 
dissipate into other forms of energy and hence contribute 
to the cosmological energy density today [19]. This con- 
tribution is called of 'vacuum realignment' . It is further 
described below. Note that the vacuum realignment con- 
tribution may be accidentally suppressed in case 1 be- 
cause the homogenized axion field happens to lie close 
to zero. 

In case 2 the axion strings radiate axions [20, 21] from 
the time of the PQ transition till t \ when the axion mass 
turns on. At t\ each string becomes the boundary of N 
domain walls. If N = 1, the network of walls bounded 
by strings is unstable [22, 23] and decays away. If N > 1 
there is a domain wall problem [24] because axion do- 
main walls end up dominating the energy density, result- 
ing in a universe very different from the one observed to- 
day. There is a way to avoid this problem by introducing 
an interaction which slightly lowers one of the N vacua 
with respect to the others. In that case, the lowest vacuum 
takes over after some time and the domain walls disap- 
pear. There is little room in parameter space for that to 
happen and we will not consider this possibility further 
here. A detailed discussion is given in Ref. [25]. Hence- 
forth, we assume N = I. 

In case 2 there are three contributions to the axion 
cosmological energy density. One contribution [20, 21, 
26, 27, 28, 29, 30] is from axions that were radiated by 
axion strings before t\. A second contribution is from 
axions that were produced in the decay of walls bounded 
by strings after /i [27, 31, 32, 25]. A third contribution is 
from vacuum realignment [19]. 

Let me briefly indicate how the vacuum alignment 



contribution is evaluated. Before time t \ , the axion field 
did not oscillate even once. Soon after t\, the axion mass 
is assumed to change sufficiently slowly that the total 
number of axions in the oscillations of the axion field 
is an adiabatic invariant. The number density of axions 
at time t \ is 

1 1 

n a {h)~-m a {h){a 2 {h))~nfl- (6) 
I t\ 

In Eq. (6), we used the fact that the axion field a{x) is ap- 
proximately homogeneous on the horizon scale t \ . Wig- 
gles in a (x) which entered the horizon long before t\ have 
been red-shifted away [33]. We also used the fact that the 
initial departure of a(x) from the nearest minimum is of 
order f a . The axions of Eq. (6) are decoupled and non- 
relativistic. Assuming that the ratio of the axion number 
density to the entropy density is constant from time t\ till 
today, one finds [19, 25] 
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for the ratio of the axion energy density to the critical 
density for closing the universe, h is the present Hubble 
rate in units of 100 km/s.Mpc. The requirement that ax- 
ions do not overdose the universe implies the constraint 
m a > 6 • 10~ 6 eV. 

The contribution from axion string decay has been de- 
bated over the years. The main issue is the energy spec- 
trum of axions radiated by axion strings. Battye and Shel- 
lard [28] have carried out computer simulations of bent 
strings (i.e. of wiggles on otherwise straight strings) and 
have concluded that the contribution from string decay 
is approximately ten times larger than that from vacuum 
realignment, implying a bound on the axion mass ap- 
proximately ten times more severe, say m a > 6 • 10~ 5 
eV instead of m a > 6 • 10 eV. My collaborators and 
I have done simulations of bent strings [27], of circular 
string loops [27, 30] and non-circular string loops [30]. 
We conclude that the string decay contribution is of the 
same order of magnitude than that from vacuum realign- 
ment. Yamaguchi, Kawasaki and Yokoyama [29] have 
done computer simulations of a network of strings in an- 
expanding universe, and concluded that the contribution 
from string decay is approximately three times that of 
vacuum realignment. The contribution from wall decay 
has been discussed in detail in ref. [25]. It is probably 
subdominant compared to the vacuum realignment and 
string decay constributions. 

It should be emphasized that there are many sources 
of uncertainty in the cosmological axion energy density 
aside from the uncertainty about the contribution from 
string decay. The axion energy density may be diluted 
by the entropy release from heavy particles which decou- 
ple before the QCD epoch but decay afterwards [34], or 



by the entropy release associated with a first order QCD 
phase transition. On the other hand, if the QCD phase 
transition is first order [35], an abrupt change of the ax- 
ion mass at the transition may increase £l a . If inflation 
occurs with reheat temperature less than Tpg, there may 
be an accidental suppression of £l a because the homog- 
enized axion field happens to lie close to a CP conserv- 
ing minimum. Because the RHS of Eq. (7) is multiplied 
in this case by a factor of order the square of the ini- 
tial vacuum misalignment angle which is randomly 
chosen between — % and +n, the probability that Q. a is 
suppressed by a factor x is of order y/x. This rule can- 
not be extended to arbitrarily small x however because 
quantum mechanical fluctuations in the axion field dur- 
ing the epoch of inflation do not allow the suppression 
to be perfect [36]. Recently, Kaplan and Zurek proposed 
a model [37] in which the axion decay constant f a is 
time-dependent, the value f a (t\ ) during the QCD phase- 
transition being much smaller than the value f a today. 
This yields a suppression of the axion cosmological en- 
ergy density by a factor {^j^-) 2 compared to the usual 
case [replace f a by f a (t\) in Eq. (6)]. 

The axions produced when the axion mass turns on 
during the QCD phase transition are cold dark matter 
(CDM) because they are non-relativistic from the mo- 
ment of their first appearance at 1 GeV temperature. 
Studies of large scale structure formation support the 
view that the dominant fraction of dark matter is CDM. 
Any form of CDM necessarily contributes to galactic ha- 
los by falling into the gravitational wells of galaxies. 
Hence, there is excellent motivation to look for axions 
as constituent particles of our galactic halo. 

Finally, let's mention that there is a particular kind of 
dumpiness [38, 25] which affects axion dark matter if 
there is no inflation after the Peccei-Quinn phase transi- 
tion. This is due to the fact that the dark matter axions are 
inhomogeneous with dp /p ~ 1 over the horizon scale at 
temperature T\ ~ 1 GeV, when they are produced at the 
start of the QCD phase-transition, combined with the fact 
that their velocities are so small that they do not erase 
these inhomogeneities by free-streaming before the time 
t eq of equality between the matter and radiation energy 
densities when matter perturbations can start to grow. 
These particular inhomogeneities in the axion dark mat- 
ter are in the non-linear regime immediately after time t eq 
and thus form clumps, called 'axion mini-clusters' [38]. 
They have mass M mc ~ 10~ 13 M© and size l mc ~ 10 13 cm. 

DARK MATTER AXION DETECTION 

An electromagnetic cavity permeated by a strong static 
magnetic field can be used to detect galactic halo axions 
[39]. The relevant coupling is given in Eq. (4). Galactic 



halo axions have velocities j3 of order 10 3 and hence 
their energies E a = m a + jm a f5 2 have a spread of order 
10~ 6 above the axion mass. When the frequency (0 = 
2nf of a cavity mode equals m a , galactic halo axions 
convert resonantly into quanta of excitation (photons) 
of that cavity mode. The power from axion — > photon 
conversion on resonance is found to be [39, 40]: 

P =(^t) VB 2 Q p a C—Mm(Q L ,Q a ) 
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where V is the volume of the cavity, Bq is the magnetic 
field strength, Q L is its loaded quality factor, Q a = 10 6 is 
the 'quality factor' of the galactic halo axion signal (i.e. 
the ratio of their energy to their energy spread), p a is the 
density of galactic halo axions on Earth, and C is a mode 
dependent form factor given by 

f v d xE a - Bo 

C = — (9) 
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where Bq(x) is the static magnetic field, E m (x)e' m is the 
oscillating electric field and e is the dielectric constant. 
Because the axion mass is only known in order of magni- 
tude at best, the cavity must be tunable and a large range 
of frequencies must be explored seeking a signal. The 
cavity can be tuned by moving a dielectric rod or metal 
post inside it. 

For a cylindrical cavity and a homogeneous longitu- 
dinal magnetic field, C = 0.69 for the lowest TM mode. 
The form factors of the other modes are much smaller. 
The resonant frequency of the lowest TM mode of a 
cylindrical cavity is /=115 MHz (^r) where R is the 
radius of the cavity. Since 10~ 6 eV = 2k (242 MHz), a 
large cylindrical cavity is convenient for searching the 
low frequency end of the range of interest. To extend 
the search to high frequencies without sacrifice in vol- 
ume, one may power-combine many identical cavities 
which fill up the available volume inside a magnet's bore 
[41, 42]. This method allows one to maintain C = 0(1) 
at high frequencies, albeit at the cost of increasing engi- 
neering complexity as the number of cavities increases. 

Axion dark matter searches were carried out at 
Brookhaven National Laboratory [43], the University of 
Florida [44], Kyoto University [45], and by the ADMX 
collaboration [46, 47, 48, 49, 50] at Lawrence Livermore 
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FIGURE 2. Axion couplings and masses excluded by the 
ADMX experiment. Also shown are the KSVZ and DFSZ 
model predictions. Indicated on the insert are the regions ex- 
cluded by the pilot experiments at Brookhaven National Labo- 
ratory (RBF) and the University of Florida (UF). All limits are 
scaled to the local halo density p a = 7.5 10~ 25 g/cm 3 [51]. 



National Laboratory. The limits published by the ADMX 
collaboration in ref. [47] are reproduced in Fig. 2. By 
definition, g arr = Updated limits can be found in 
refs. [48,49]. 

The ADMX experiment is equipped with a high reso- 
lution spectrometer which allows us to look for narrow 
peaks in the spectrum of microwave photons caused by 
discrete flows, or streams, of dark matter axions in our 
neighborhood. In many discussions of cold dark matter 
detection it is assumed that the distribution of CDM par- 
ticles in galactic halos is isothermal. However, there are 
excellent reasons to believe that a large fraction of the 
local density of cold dark matter particles is in discrete 
flows with definite velocities [52]. Indeed, because CDM 
has very low primordial velocity dispersion and negligi- 
ble interactions other than gravity, the CDM particles lie 
on a 3-dim. hypersurface in 6-dim. phase-space. This im- 
plies that the velocity spectrum of CDM particles at any 
physical location is discrete, i.e., it is the sum of distinct 
flows each with its own density and velocity. 

We searched for the peaks in the spectrum of mi- 
crowave photons from axion to photon conversion that 
such discrete flows would cause in the ADMX detector. 
We found none and placed limits [50] on the density of 
any local flow of axions as a function of the flow ve- 
locity dispersion over the axion mass range 1.98 to 2.17 
jtieV. Our limit on the density of discrete flows is approx- 
imately a factor three more severe than our limit on the 
total local axion dark matter density. 

The ADMX experiment is presently being upgraded 
to replace the HEMT (high electron mobility transis- 
tors) receivers we have used so far with SQUID mi- 
crowave amplifiers. HEMT receivers have noise temper- 
ature T n ~ 3 K [53] whereas T n <~ 0.05 K was achieved 
with SQUIDs [54]. In a second phase of the upgrade, 



the experiment will be equipped with a dilution refriger- 
ator to take full advantage of the lowered electronic noise 
temperature. When both phases of the upgrade are com- 
pleted, the ADMX detector will have sufficient sensitiv- 
ity to detect DFSZ axions at even a fraction of the local 
halo density. 



SOLAR AXION DETECTION 

The conversion of axions to photons in a magnetic 
field can also be used to look for solar axions [39, 
55, 56, 57, 58]. The flux of solar axions on Earth is 
L±i2 T-(7rfe) 2 ( i %) 2 from the Primakoff conversion of 

sec cm 2 v 0.36 / v eV ' 

thermal photons in the Sun [55]. The actual flux may 
be larger because other processes, such as Compton-like 
scattering, contribute if the axion has an appreciable cou- 
pling to the electron. At any rate the flux is huge com- 
pared to what can be produced by man-made processes 
on Earth and it is cost free. Solar axions have a broad 
spectrum of energies of order the temperature in the so- 
lar core, from one to a few keV. 

Since the magnetic field is homogeneous on the length 
scale set by the axion de Broglie wavelength, the final 
photon is colinear with the initial axion. The photon and 
axion also have the same energy assuming the magnetic 
field is time-independent. The a — > y conversion proba- 
bility is [39, 55] 



? arY ) 2 (B 0± LF(q)) 2 



(10) 



if Bq_i_(z) = Boj_b(z) is the magnetic field transverse to 
the direction of the colinear axion and photon, z is the 
coordinate along this direction, L is the depth over which 
the magnetic field extends and F(q) is the form factor 

F (l) = j\ j L dze' qz b(z) | (11) 
L Jo 
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where q = k y -k a =E a - ^/E 2 -m 2 
mentum transfer. If the magnetic field is homogeneous 
(b=l), then 
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For qL ^> 1, the conversion probability goes as sin 2 (-^) 
because the axion and photons oscillate into each other 
back and forth. The form factor F (q) can be improved by 
filling the conversion region with a gas whose pressure is 
adjusted in such a way that the plasma frequency, which 
acts as an effective mass for the photon, equals the axion 
mass [55]. 
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FIGURE 3. Regions of (m„ , g a yy) space ruled out by various 
solar axion searches. 



Multiplying the flux times the conversion probability, 
one obtains the event rate: 
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The final state photons are soft x-rays which may be de- 
tected with good efficiency. There are radioactive back- 
grounds to worry about, however. 

The above type of detector is usually referred to as 
an axion helioscope. If a signal is found due to axions 
or familons, the detector becomes a marvelous tool for 
the study of the solar interior. Experiments were carried 
out at Brookhaven National Lab. [56], the University of 
Tokyo [57], and more recently at CERN by the CAST 
collaboration [58]. Fig. 3 shows the recently published 
CAST limits, as well as previous limits from the Tokyo 
group, and from the DAMA [59] and SOLAX [60] col- 
laborations. The latter two experiments searched for so- 
lar axions by looking for their Primakoff conversion to 
photons in Nal and Ge crystals respectively. 



LASER EXPERIMENTS 



magnet, a few of the photons convert to axions. Another 
dipole magnet is set up in line with the first, behind a 
wall. Since the axions go through the wall unimpeded, 
this setup allows one to 'shine light through the wall.' 
An experiment of this type was carried out by the RBF 
collaboration [62]. Also a proposal has been made to 
use decommissioned HERA dipole magnets [63] for this 
purpose. Compared with a solar axion search, the 'shin- 
ing light through walls' experiment has the advantage 
of greater control over experimental parameters. But the 
signal is much smaller because one pays twice the price 
of the very small axion-photon conversion rate. 

It was also proposed [64] to look for the effect of the 
axion on the propagation of light through a magnetic 
field. If the photon beam is linearly polarized and the po- 
larization direction is at an angle to the direction of the 
magnetic field, the plane of polarization turns because 
the polarization component parallel to the magnetic field 
gets depleted whereas the perpendicular component does 
not. There is an additional effect of birefringence because 
the component of light polarized parallel to the magnetic 
field mixes with the axion and hence moves more slowly 
than in vacuo. Birefringence affects the ellipticity of the 
polarization as the light travels on. Experiments search- 
ing for polarization rotation and birefringence of laser 
light passing through a magnetic field were carried out 
by the the RBF [65] and PVLAS [66] collaborations. In 
both experiments, the optical path was boosted by pass- 
ing the laser beam many times through an optical cav- 
ity within the magnet. A positive result has very recently 
been reported by the PVLAS collaboration. 

The PVLAS collaboration [66] finds an optical rota- 
tion of (3.9 ±0.5) 10~ 12 radians per pass, at 5 Tesla. 
When interpreted in terms of an axion-like particle a with 
coupling to two photons r 5£ a yy — jf a E ' B, the range of 
mass and coupling consistent with the PVLAS effect is 
0.7 meV < m a < 2 meV and 1 • 10 5 GeV <M a <6- 10 5 
GeV At face value the PVLAS result disagres with the 
limits on axion-like particles from stellar evolution and 
solar axion searches, since -X- is the same as g u yy in 
Fig. 3. On the other hand, there may be ways to make the 
various results consistent with one another, as proposed 
for instance in ref . [67] . 



Eqs. (10,12) give the conversion probability in a static 
magnetic field of an axion to a photon of the same en- 
ergy. The polarization of the photon is parallel to the 
component of the magnetic field transverse to the di- 
rection of motion. The inverse process, conversion of 
such a photon to an axion, occurs with the same prob- 
ability p. K. van Bibber et al. [61] proposed a 'shining 
light through walls' experiment in which a laser beam 
is passed through a long dipole magnet like those used 
for high-energy physics accelerators. In the field of the 
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